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We critically examine the gauge, and field-parametrization dependence of renormalization group 
flows in the vicinity of non-Gaufiian fixed points in quantum gravity. While physical observables 
are independent of such calculational specifications, the construction of quantum gravity field the¬ 
ories typically relies on off-shell quantities such as jj functions and generating functionals and thus 
face potential stability issues with regard to such generalized parametrizations. We analyze a two- 
parameter class of covariant gauge conditions, the role of momentum-dependent field rescalings and 
a class of field parametrizations. Using the product of Newton and cosmological constant as an 
indicator, the principle of minimum sensitivity identifies stationary points in this parametrization 
space which show a remarkable insensitivity to the parametrization. In the most insensitive cases, 
the quantized gravity system exhibits a non-Gaufiian UV stable fixed point, lending further support 
to asymptotically free quantum gravity. One of the stationary points facilitates an analytical deter¬ 
mination of the quantum gravity phase diagram and features ultraviolet and infrared complete RG 
trajectories with a classical regime. 


I. INTRODUCTION 

Physical observables are independent of their compu¬ 
tational derivation. Still, many practical computations 
are based on convenient choices for intermediate auxil¬ 
iary tools such as coordinate systems, gauges, etc. Ap¬ 
propriate parametrizations of the details of a system sim¬ 
ply decrease the computational effort. Beyond pure effi¬ 
ciency aspects, such suitable parametrizations can also be 
conceptually advantageous or even offer physical insight. 
This is similar to coordinate choices in classical mechan¬ 
ics where polar coordinates with respect to the ecliptic 
plane in celestial mechanics support a better understand¬ 
ing in comparison with, say, Cartesian coordinates with 
a z axis pointing towards Betelgeuse. 

Appropriate parameterizations become particularly 
significant in quantum calculations. While on-shell 
quantities such as 5'-matrix elements are invariant ob¬ 
servables [1-3], off-shell quantities generically feature 
parametrization dependencies, including gauge-, field- 
parametrization and regularization-scheme dependencies 
[4-6]. Further ordering schemes such as perturbative 
expansions may defer such dependencies to higher or¬ 
ders (such as scheme dependence in mass-independent 
schemes), but these are merely special and not always 
useful limits. Approximation schemes that can also deal 
with non-perturbative regimes may even introduce fur¬ 
ther artificial parametrization dependencies which have 
to be carefully removed (e.g., discretization artefacts in 
lattice regularizations). 

In an ideal situation, this parametrization dependence 
of a nonperturbative approximation could be quanti¬ 
fied and proven to be smaller than the error of the 
truncated solution. However, as soon as a result is 
parametrization dependent, it is likely that some patho¬ 
logical parametrization can be constructed that modifies 
the result in an arbitrary fashion. This suggests to look 
for general criteria of good parametrizations that mini¬ 
mize the artificial dependence in approximation schemes 
which adequately capture the physical mechanisms. 


A-priori criteria suggest the construction of 
parametrizations that support the identification of 
physically relevant degrees of freedom, such as the use 
of Coulomb-Weyl gauge in quantum optics, or the use 
of pole-mass regularization schemes in heavy-quark 
physics. Further a-priori criteria include symmetry 
preserving properties (covariant gauges, non-linear 
field parametrizations) or strict implementations of a 
parametrization condition such as the Landau-gauge 
limit a —>■ 0. A major advantage of the latter is that 
some redundant degrees of freedom decouple fully from 
the dynamical equations in such a limit. 

Good parametrizations may also be identified a poste¬ 
riori by allowing for a family of parametrizations and 
identifying stationary points in the parameter space. 
This realizes the principle of minimum sensitivity [7, 8] 
(originally advocated for regularization-scheme depen¬ 
dencies), suggesting those points as candidate parame¬ 
ters for minimizing the influence of parametrization de¬ 
pendencies. 

In the present work, we investigate a two-parameter 
family of covariant gauges, a family of field parametriza¬ 
tions and the role of momentum-dependent field rescal¬ 
ings in quantum gravity in this spirit. The family of 
gauges includes a (non-harmonic) generalization of the 
harmonic gauge (De-Donder gauge), the latter being par¬ 
ticularly useful for the analysis of gravitational waves 
which presumably are the asymptotic states of quantum 
gravity. The a-priori criteria suggest to implement this 
gauge in the Landau-gauge limit to decouple a redundant 
part of the Hilbert space. In fact, in this limit we find a 
subtlety in the form of a degeneracy in the subspace of 
scalar field components which is special to gravity. 

We also investigate a one-parameter family of field 
parametrizations that includes the most widely used lin¬ 
ear split [9] as well as the exponential split [10-14] studied 
more recently in the context of asymptotic safety [15-18] 
- both of which find support by discriminative a-priori 
arguments. We also take a brief look at the most general 
ultra-local four-parameter family of parametrizations to 
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quadratic order, corroborating the results of the one- 
parameter family. In addition, we study the influence 
of momentum-dependent field rescalings which are com¬ 
monly used in gravity in connection with the York decom¬ 
position. In the context of the functional renormalization 
group (FRG) [19] which provides a tool to study quantum 
gravity nonperturbatively [20] , these parametrization de¬ 
pendencies can mix nontrivially with the regularization 
of the spectrum of fluctuations. Therefore, the analysis 
of parametrization dependencies also explores implicitly 
the stability of the system in the ultraviolet. 

Interestingly, we observe a nontrivial interplay between 
all these parametrization dependencies. Still, several sta¬ 
tionary points can be observed in the results for the RG 
flow where the system develops a remarkable insensitiv¬ 
ity to the details of the parametrization choices. In par¬ 
ticular, for the stable parametrizations, we observe the 
existence of a UV stable non-Gaufiian fixed point which 
provides further quantitative evidence for the existence of 
an asymptotically safe metric quantum gravity [21, 22] . 
In the stationary regime of the parametrization based 
on the exponential split, the resulting RG flow exhibits 
several remarkable properties: (1) a possible dependence 
on the residual gauge parameter drops out implying an 
enhanced degree of gauge invariance, (2) the RG flow be¬ 
comes particularly simple, such that the phase diagram 
in the plane of Newton and cosmological constant can 
be computed analytically, (3) no singularities arise in the 
flow, such that a large class of RG trajectories (including 
those with a classical regime) can be extended to arbi¬ 
trarily high and low scales, (4) the UV critical exponents 
are real and close to their canonical counterparts, and (5) 
indications are found that the asymptotic safety scenario 
may not extend straightforwardly to dimensions much 
higher than d = 4. 


II. QUANTUM GRAVITY AND 
PARAMETRIZATIONS 

The technical goal of quantum gravity is to construct 
a functional integral over suitable integration variables 
which in the long-range limit can be described by a 
diffeomorphism-invariant effective field theory of metric 
variables approaching a classical regime for a wide range 
of macroscopic scales. The fact that the first part of this 
statement is rather unspecific is reflected by the large 
number of legitimate quantization proposals [23-25]. In¬ 
dependently of the precise choice of integration variables, 
a renormalization group approach appears useful in order 
to facilitate a scale-dependent description of the system 
and a matching to the long-range classical limit which is 
given at least to a good approximation by an (effective) 
action of Einstein-Hilbert type: 


Here, we have already introduced a momentum scale k, 
expressing the fact that this effective description should 
a priori hold only for a certain range of classical scales. 
In this regime, we have Zji = l/(167rGN) with the New¬ 
ton constant Gn, and A parametrizing the cosmological 
constant. In a quantum setting, Zn plays the role of 
a (dimensionful) wave-function renormalization, and Gn 
and A are expected to be replaced by their running coun¬ 
terparts depending on the scale k. 

In the present work, we confine ourselves to a quantum 
gravity field theory assuming that the metric itself is al¬ 
ready a suitable integration variable. A first step towards 
a diffeomorphism-invariant functional integral then pro¬ 
ceeds via the Faddeev-Popov method involving a gauge 
choice for intermediate steps of the calculation. In this 
work, we use the background-field gauge with the gauge¬ 
fixing quantity, 

( 2 ) 

which should vanish if the gauge condition is exactly 
matched. Here, ga/s is the full (fluctuating) met¬ 
ric, whereas ga^ denotes a fiducial background metric 
which remains unspecified, but assists to keep track of 
diffeomorphism invariance within the background-field 
method. Gauge-fixing is implemented in the functional 
integral by means of the gauge-fixing action 

Tgf = § / (3) 

More precisely, this gauge choice defines a two-parameter 
(a, /3) family of covariant gauges. For instance, the choice 
/I = 1 corresponds to the harmonic/De-Donder gauge 
which together with a = 1 (Feynman gauge) yields a va¬ 
riety of technical simplifications, being used in standard 
effective field theory calculations [26-28] as well as in 
functional RG studies [20, 29] of quantum gravity. More 
conceptually, the Landau-gauge limit a —>■ 0 appears fa¬ 
vorable, as it implements the gauge condition in a strict 
fashion and thus should be a fixed point under RG evo¬ 
lution [30, 31] . 

In the Euclidean formulation considered here, the pa¬ 
rameter a is bound to be non-negative to ensure the pos¬ 
itivity of the gauge-fixing part of the action (this restric¬ 
tion may not be necessary for a Lorentzian formulation). 
The parameter j3 can be chosen arbitrarily except for the 
singular value /3sing = d — 1. To elucidate this singularity, 
let us take a closer look at the induced Faddeev-Popov 
ghost term: 

Pgh = - y d‘'a:VIG^^f^G^ = —, (4) 

where v‘' characterizes the vector field along which we 
study the Lie derivative generating the coordinate trans¬ 
formations, 

dgap 


^k = -J d^xV5Zi^(i?-2A). 


( 1 ) 


Sv‘ 


S „ 6 

-^Zygap - 2 — D(^aVp)- 


( 5 ) 
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The corresponding variation of the gauge-fixing condition 
yields 

Let us decompose the vector Svjs into a transversal part 
6vJ and a longitudinal part Dp5x- For the following 
argument, it suffices to study the limit of the quantum 
metric approaching the background metric > gap, 

which diagrammatically corresponds to studying the in¬ 
verse ghost propagator ignoring higher vertices 

+ ^((d - l-p)D^^D,+AR>:;)D''5x + 0{g - g). 

In this form it is obvious that the longitudinal direction 
D'^Sx is not affected by the gauge fixing for /3 = d — 1 
to zeroth order in the curvature. In other words, the 
gauge fixing is not complete for this singular case /3sing = 
d — 1. This singularity is correspondingly reflected by 
the ghost propagator. The Faddeev-Popov operator in 
Eq. (4) reads 

= 2r^D<^D^^gp), - 2^-^r^D^^Dagp,. ( 8 ) 

Decomposing the ghost fields (7^, also into transversal 
Cj, and longitudinal parts D^fj, D''g we find for the 
ghost Lagrangian, 

= Cj + Ri^) (9) 

-n D^FR^''D^D}jg + 0{g - g). 

where we have performed partial integrations in order to 
arrive at a convenient form and dropped covariant deriva¬ 
tives of the curvature. This form of the inverse propaga¬ 
tor of the ghosts makes it obvious that a divergence of 
the form arises in the longitudinal parts. This di¬ 

vergence at /3sing = d — 1 related to an incomplete gauge 
fixing will be visible in all our results below. 

Let us now turn to the metric modes. As a techni¬ 
cal tool, we parametrize the fully dynamical metric g^v 
in terms of a fiducial background metric gf^, and fluc¬ 
tuations h^i, about the background. Background inde¬ 
pendence is obtained by keeping g^^, arbitrary and re¬ 
quiring that physical quantities such as scattering ampli¬ 
tudes are independent of Still, these requirements 
do not completely fix the parametrization of the dynam¬ 
ical field g = g[g] h]. Several parametrizations have been 
used in concrete calculations. The most commonly used 
parametrization is the linear split [9] 

gfiu = gfiu + hfj,,,. ( 10 ) 

By contrast, the exponential split [10-14] 

gi.v =gt,p{e’^Yu, (11) 


is a parametrization that has been discussed more re¬ 
cently to a greater extent [15-18]. In both cases, h is 
considered to be a symmetric matrix field (with indices 
raised and lowered by the background metric). If a path 
integral of quantum gravity is now defined by some suit¬ 
able measure T>h, it is natural to expect that the space 
of dynamical metrics g is sampled differently by the two 
parametrizations, implying different predictions at least 
for off-shell quantities - unless the variable change from 
(10) to (11) is taken care of by suitable (ultralocal) Ja- 
cobians. While a parametrization (and gauge-condition) 
independent construction of the path integral has been 
formulated in a geometric setting [9, 32-35], its usability 
is hampered by the problem of constructing the full de¬ 
composition of h in terms of fluctuations between physi¬ 
cally inequivalent configurations and fluctuations along 
the gauge orbit. Geometric functional RG flows have 
been conceptually developed in [36], with first results 
for asymptotic safety obtained in [37], and recently to 
a leading-order linear-geometric approximation in [38]. 
The relation between the geometric approach and the 
exponential parametrization was discussed in [16]. 

In the present work, we take a more pragmatic view¬ 
point, and consider the different parameterizations of 
Eqs. (10) and (11) as two different approximations of an 
ideal parametrization. Since the functional RG actually 
requires the explicit form of g[g', h] only to second order 
in h (in the single-metric approximation, see below), we 
mainly consider a one-parameter class of parametriza¬ 
tions of the type 

gpi/ = g^L/ + + 2 0{h^). ( 12 ) 

For T = 0, we obtain the linear split, whereas r = 1 is 
exactly related to the exponential split within our trun¬ 
cation. Incidentally, it is straightforward to write down 
the most general, ultra-local parametrization to second 
order that does not introduce a scale, 

g^u — g^v T 

+ 0{hY. (13) 

Here, h = hf. is the trace of the fluctuation. As men- 
tioned above, third and higher-order terms will not con¬ 
tribute to our present study anyway. Instead of exploring 
the full parameter dependence, we will highlight some in¬ 
teresting results in this more general framework below. 

The key ingredient for a quantum computation is the 
propagator of the dynamical field. In our setting, its 
inverse is given by the second functional derivative (Hes¬ 
sian) of the action (1) including the gauge fixing (3) with 
respect to the fluctuating field h, 
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1 p(2)Kl/ 

Zr h= 0 ,C =0 

= - [8a - (1 + /3)^]r"5a/5) (-C^) 

- 

+ ^ ^ 4 ~ ir'^Df^cDp) + gc.pD^^D''^) 

+ i(2(l-r)<5S^-r"M(i?-2Afc) 

- (1 - r)i?|:<5;j + + i?a/9r") - R^cTp), 

(14) 


Here and in the following, we specialize to d = 4, except 
if stated otherwise. A standard choice for the gauge pa¬ 
rameters is harmonic DeDonder gauge with a = 1 = /3 
for which the second and third lines simplify consider¬ 
ably. Simplifications also arise for the exponential split 
r = 1; in particular, a dependence on the cosmological 
constant remains only in the trace mode ~ g'^'^ga/ 3 - 
A standard tool for dealing with the tensor structure 
of the propagator is the York decomposition of the fluc¬ 
tuations h^i, into transverse traceless tensor modes, a 
transverse vector mode and two scalar modes, 

hfj,u = h— —gfiuD ^ '^g^tyh, 

(15) 

D^hl, = 0, r''hl, = Q, = (16) 

It is convenient to split into a pure kinetic part V 
which has a nontrivial flat-space limit, and a curvature- 
dependent remainder = 0{R). This facilitates an 
expansion of the propagator (F^^^)”^ = {V -F)“^ = 

OO 

£ (-ip-ijc-)"ip-i. 

n—0 

Let us first concentrate on the kinetic part V'. 


Whereas the choice of a and /I, in principle, are inde¬ 
pendent, there can arise a subtle interplay with certain 
regularization strategies as will be highlighted in the fol¬ 
lowing. 

By taking the limit a —>■ 0 while keeping /3 finite, we 
make the gauge fixing explicit, especially we And for the 
gauge-dependent modes 




-.-1 


5^ 

'-'a. 1 


_ wa-2 


(ah)^ (^A-( 3-/32 + ( 3 +/ 32 )r)Afc 


( 20 ) 

' /32 -6/3A\ 

-6^A 36A2 ) ■ 

( 21 ) 


The transverse mode decouples linearly with a —>■ 0 
and hence is pure gauge in the present setting. Whereas 
finite parts seem to remain in the (cr/i) subspace, we ob¬ 
serve that the matrix in (21) becomes degenerate in 
this limit (e.g., the determinant of the matrix in Eq. (21) 
is zero). Effectively, only one scalar mode remains in the 
propagator. The nature of this scalar mode is a function 
of the second gauge parameter: taking the limit /3 —>■ oo, 
the remaining scalar mode can be identified with tr, while 
the limit (5^0 leaves us with a pure h mode. 

Whereas the transverse modes in Eq. (20) decouple 
smoothly in the limit a —t 0, the decoupling of the scalar 
mode in Eq. (21) is somewhat hidden in the degeneracy 
of the scalar sector with the corresponding eigenmode 
depending on /3. This can lead to a subtle interplay 
with regularization techniques for loop diagrams as can 
be seen on rather general grounds by the following argu¬ 
ment. Structurally, the propagator in the (u/i) sector has 
the following form in the limit a —>■ 0 and for small but 
finite P, cf. Eq. (21) 

(V, A-i ^ r22l 

[R(<7h)) ^ 0(1) j- 

Regularizations of traces over loops built from this prop¬ 
agator are typically adjusted to the spectrum of the in¬ 
volved operators. Let us formally write this as 


^.T-„, = ^<5-(A-2(l-r)A.), 

^5T'l = ^<5((A(A-2a(l-r)Afc), 


(17) 

(18) 


_ /3. ( 3 -<-)A-4 £ ( i ^ A y A^ 

(^/«) A.{P-a)A^ (/3--3a)A+4.(l+.)A. ^ 


(19) 


where A = —D^. In this form it is straightforward 
to calculate the propagator {V)~^. In particular, the 
transverse traceless mode does not exhibit any de¬ 
pendence on the gauge parameters. As discussed in the 
introduction, a-priori criteria suggest the Landau-gauge 
limit a —>■ 0 as a preferred choice for the gauge fixing, 
as it strictly implements the gauge-fixing condition. It 
thus should also be a fixed point of the RG flow [30, 31]. 


Tr 


CnV-\...) 


(23) 


where Cr denotes a regularizing operator and the ellip¬ 
sis stands for further vertices and propagators. Now, it is 
often useful to regularize all fluctuation operators at the 
same scale, e.g., the spectrum of all A’s should be cut off 
at one and the same scale kP. Therefore, the regularizing 
operator C-r inherits its tensor structure from the Hes¬ 
sian F(2) of Eq. (14). In the ( ah) sector, the regularizing 
operator can hence acquire the same dependence on the 
gauge-parameters as in Eq. (19), 


^TZ,{<7h) 


I foil) OiP)\ 

a 'v(^(^) OiP^))^ 


(24) 


for a —>■ 0 and small /?. The complete scalar contribution 
to traces of the type (23) would then be of the parametric 
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form, 


W). (25) 

For finite /3, such regularized traces can thus be afflicted 
with divergencies in the Landau-gauge limit a —> 0. If 
this happens, we still have the option to choose suitable 
values of /?. In fact, Eq. (25) suggest that still a whole 
one-parameter family of gauges exists in the Landau- 
gauge limit, if we set /3 = 7 • -^/a, with arbitrary real but 
finite gauge parameter 7 distinguishing different gauges. 

We emphasize that this is a rather qualitative analysis. 
Since the limit of products is not necessarily equal to the 
product of limits, the trace over the matrix structure 
of the above operator products can still eliminate this 
l/cK divergence, such that any finite value of /3 remains 
admissible. 

In the following we observe that the appearance of the 
1 ja. divergence depends on the explicit choice of the regu¬ 
larization procedure, as expected. Still, as this discussion 
shows, even if this divergence occurs, it can perfectly well 
be dealt with by choosing /3 = and still retaining 

a whole one-parameter family of gauges in the Landau 
gauge limit. 


III. GRAVITATIONAL RG FLOW 

For our study of generalized parametrization depen¬ 
dencies of gravitational RG flows, we use the func¬ 
tional RG in terms of a flow equation for the effective 
average action (Wetterich equation) [19] amended by 
the background-field method [39-41] and formulated for 
gravity [ 20 ] 


dtrk[g,g] = ^STr +7^fe) ' 


dt = k- 


'dk' 
(26) 

Equation (26) describes the flow of an action functional 
Ffc as a function of an RG scale k that serves as a regu- 
larization scale for the infrared fluctuations. Here, F ),' 
denotes the Hessian of the action with respect to the fluc¬ 
tuation field g, at fixed background g. The details of the 
regularization are encoded in the choice of the regulator 
TZk- Suitable choices of TZk guarantee that F^ becomes 
identical to the full quantum effective action in the limit 
A: —>■ 0, and approaches the bare action for large scales 
k —)■ Auv 00 (where Auv denotes a UV cutoff). For 
reviews in the present context, see [29, 42-47]. 

Whereas exact solutions of the flow equation so far 
have only been found for simple models, approximate 
nonperturbative flows can be constructed with the help 
of systematic expansion schemes. In the case of gravity, 
a useful scheme is given by expanding F^ in powers of 
curvature invariants. The technical difficulties then lie in 
the construction of the inverse of the regularized Hessian 


+ TZk) corresponding to the regularized propaga¬ 

tor, and performing the corresponding traces (the super¬ 


trace STr includes a minus sign for Grassmann degrees 
of freedom, i.e., Faddeev-Popov ghosts). 

A conceptual difficulty lies in the fact that rk[g,g] 
should be computed on a subspace of action functionals 
that satisfy the constraints imposed by diffeomorphism 
invariance and background independence. In general, 
this requires to work with g and g independently during 
large parts of the computation [48-51]. Such bi-metric 
approaches can, for instance, be organized in the form 
of a vertex expansion on a flat space as put forward re¬ 
cently in [52-54], or via a level expansion as developed 
in [55], see [56-58] for further bi-metric results. For the 
present study of parametrization dependencies, we con¬ 
fine ourselves to a single-metric approximation, defined 
by identifying g with g on both sides of the flow equa¬ 
tion, after the Hessian has been analytically determined. 
In the following, we therefore do no longer have to distin¬ 
guish between the background field and the fluctuation 
field as far as the presentation is concerned, and hence 
drop the bar notation for simplicity. 

Spanning the action in terms of the Einstein-Hilbert 
truncation ( 1 ) and neglecting the flow of the gauge-fixing 
and ghost sector [59-61], we use the universal RG ma¬ 
chine [62-64] as our computational strategy. The key 
idea is to subdivide the Hessian F), into a kinetic part 
and curvature parts with a subsequent expansion in the 
curvature. This is complicated by terms containing un- 
contracted covariant derivatives in F),' which could in¬ 
validate the counting scheme. Within the present trun¬ 
cation, this problem is solved with the aid of the York 
decomposition (15). This helps both to set up the cur¬ 
vature expansion as well as to invert the kinetic terms in 
the corresponding subspaces of TT, T and scalar modes. 
From a technical point of view, we use the package xAct 
[65-70] to handle the extensive tensor calculus. 

Schematically, the flow equation for the Einstein- 
Hilbert truncation can then be written as 

dtTk = j d^xv^ ( 5 TT + 5^ + 5'^'“ + + 5 J""), 

(27) 

where the first three terms denote the contributions from 
the graviton fluctuations as parametrized by the York de¬ 
composition (15). The fourth term 5®^ arises from the 
Faddeev-Popov ghost fluctuations, cf. Eq. (10). The 
last term comes from the use of transverse decom¬ 
positions of the metric (15) and the ghost fields (10). 
The corresponding functional integral measure over the 
new degrees of freedom involves Jacobians which - upon 
analogous regularization - contribute to the flow of the 
effective average action. 

At this point, we actually have a choice that serves as 
another source of parametrization dependencies studied 
in this work: one option is to formulate the regularized 
path integral in terms of the decomposed fields as intro¬ 
duced above. In that case, the Jacobians are nontrivial 
and their contribution is listed in Eq. (A21). Alter¬ 
natively, we can reintroduce canonically normalized fields 
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by means of a nonlocal field redefinition [71, 72], 



VA-Ric 

(28) 

a ^ a, 

(29) 

VA ?7 —t ?7, 

(30) 


and analogously for the longitudinal anti-ghost field fj. 
(Here, we have used (Ric^)'^ = This field re¬ 

definition goes along with another set of Jacobians con¬ 
tributing to the measure of the rescaled fields. As shown 
in [72], the Jacobians for the original York decomposition 
and the Jacobians from the field redefinition (30) cancel 
at least on maximally symmetric backgrounds. The lat¬ 
ter choice of backgrounds is sufficient for identifying the 
flows in the Einstein-Hilbert truncation. Therefore, if we 
set up the flow in terms of the redefined fields (30), the 
last term in Eq. (27) vanishes, = 0. 

For an exact solution of the flow, it would not mat¬ 
ter whether or not a field redefinition of the type (30) is 
performed. Corresponding changes in the full propaga¬ 
tors would be compensated for by the (dis-) appearance 
of the Jacobians. For the present case of a truncated 
nonperturbative flow, a dependence on the precise choice 
will, however, remain, which is another example for a 
parametrization dependence. This dependence also arises 
from the details of the regularization. The universal RG 
machine suggests to construct a regulator TZk such that 
the Laplacians A appearing in the kinetic parts are re¬ 
placed by 


A^A-bRfe(A), (31) 


where Rk{x) is a (scalar) regulator function that pro¬ 
vides a finite mass-like regularization for the long-range 
modes, e.g., Rk{x) —)■ for x but leaves the 

UV modes unaffected, Rk{x) —t 0 for a: Since 

the field redefinition (30) is nonlocal, it also affects the 
kinetic terms and thus takes influence on the precise man¬ 
ner of how modes are regularized via Eq. (31). In other 
words, the dependence of our final results on using or 
not using the field redefinition (30) is an indirect probe 
of the regularization-scheme dependence and thus of the 
generalized parametrization dependence we are most in¬ 
terested in here. 

In this work, we focus on the RG flow of the effec¬ 
tive average action parametrized by the operators of the 
Einstein-Hilbert truncation (1). For this, we introduce 
the dimensionless versions of the gravitational coupling 
and the cosmological constant. 


lenZji 


k^G, 



(32) 


and determine the corresponding RG /3 functions for g 
and A, by computing the S terms on the right hand side 
of the flow (27) to order R in the curvature. Many higher- 
order computations have been performed by now [73-84], 
essentially confirming and establishing the simple picture 
visible in the Einstein-Hilbert truncation. 


We are particularly interested in the existence of fixed 
points 5 * and A* of the jS functions, defined by 

dtg = g = I3g{g*,x*) = 0, dtX = X = l3\{g*,Xt) = 0. 

(33) 

In addition to the Gaufiian fixed point g* = 0 = A*, 
we search for a non-Gaufiian interacting fixed point, the 
existence of which is a prerequisite for the asymptotic- 
safety scenario. Physically viable fixed points should 
have a positive value for the Newton coupling and should 
be connectable by an RG trajectory with the classical 
regime, where the dimensionful couplings are approxi¬ 
mately constant, i.e., the dimensionless versions should 
scale as g ~ A ~ 1/k^. The asymptotic-safety sce¬ 
nario also requires that a possible non-Gaufiian fixed 
point has finitely many UV attractive directions. This 
is quantified by the number of positive critical exponents 
di which are defined as (— 1 ) times the eigenvalues of the 
stability matrix 9/3(g ^)/9(g, A). 

Whereas the fixed-point values g* and A* are RG 
scheme-dependent, the critical exponents 9^ are univer¬ 
sal and thus should be parametrization independent in 
an exact calculation. Also, the product 5 *A* has been 
argued to be physically observable in principle and thus 
should be universal [72] . Testing the parametrization de¬ 
pendence of the critical exponents 9i and g^A* therefore 
provides us with a quantitative criterion for the reliability 
of approximative results. 


IV. GENERALIZED PARAMETRIZATION 
DEPENDENCE 

With these prerequisites, we now explore the 
parametrization dependencies of the following scenarios: 
we consider the linear ( 10 ) and the exponential ( 11 ) split, 
both with and without field redefinition (30), and study 
the corresponding dependencies on the gauge parameters, 
focusing on a strict implementation of the gauge-fixing 
condition a —>■ 0 (Landau gauge). As suggested by the 
principle of minimum sensitivity, we look for stationary 
points as a function of the remaining parameter(s) where 
universal results become most insensitive to these gener¬ 
alized parametrizations. For the following quantitative 
studies, we exclusively use the piecewise linear regulator 
[85, 86 ], Rk{x) = (fc^ — x^)9{k'^ — x^), for reasons of sim¬ 
plicity. Studies of regulator-scheme dependencies which 
can also quantify parametrization dependencies have first 
been performed, e.g., in [72, 78] . 


A. Linear split without field redefinition 

Let us start with the case of the linear split (10) with¬ 
out field redefinition (30). Here, the degeneracy in the 
sector of scalar modes interferes with the regularization 
scheme, as illustrated in Eq. (25). Hence, in the Landau- 
gauge limit a —> 0 , we choose /3 = 7 y/a, which removes 
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Figure 1. Linear split without field redefinition: residual 
dependence of our estimates for the universal quantities on 
gauge parameter 7 in the limit a —>■ 0. We find a common 
stationary point at 7 = 0 and a remarkably small variation of 
the results on the level of 0 . 1 % for p,A* and 1 . 6 % for ReS in 
the range 7 € [— 2 , 2 ]. 


any artificial divergence, but keeps 7 as a real parameter 
that allows for a quantification of remaining parametriza- 
tion/gauge dependence. We indeed find a non-Gaufiian 
fixed point for wide range of values of 7 . The 

critical exponents form a complex conjugate pair. The 
estimates for the universal quantities and the real 
part of the 9’s (being the measure for the RG relevance 
of perturbations about the fixed point) are depicted in 
Fig. 1. 

We observe a common point of minimum sensitivity 
at 7 = 0. In a rather wide range of gauge parameter 
values 7 G [—2,2], our estimates for g*A» and Re61 vary 
only very mildy on the level of 0.1% and 1.6%. Given the 
limitations of the present simple approximation, this is a 
surprising degree of gauge independence lending further 
support to the asymptotic-safety scenario. The extrem- 
izing values at 7 = 0 are near the results of [62, 76, 77[ 
where the same gauge choice (a = /3 = 0) was used. The 
main difference can be traced back to the fact that our 
inclusion of the (dimensionful) wave-function renormal¬ 
ization in the gauge fixing term (3) renders the gauge 
parameter a dimensionless as is conventional. If we ig¬ 
nored the resulting dimensional scaling, our extremizing 


parametrization 

9* 

A. 

g*X, 

e 

nfr T — a — 'y — 0 

0.879 

0.179 

0.157 

1.986 ± i 3.064 

nfr T = 0, a = /3 = l 

0.718 

0.165 

0.119 

1.802 ± i 2.352 

frr = Q: = 0, /?=! 

0.893 

0.164 

0.147 

2.034 ± i 2.691 

frr = 0, q: = /3 = 1 

0.701 

0.172 

0.120 

1.689 ± i 2.486 

frr = Q: = 0, /? = oo 

0.983 

0.151 

0.148 

2.245 ± i 2.794 

II 

II 

8 

3.120 

0.331 

1.033 

4, 2.148 

fr r = 1.22, a = 0, /3 = 00 

3.873 

0.389 

1.508 

3.957, 1.898 


Table I. Non-Gaufiian fixed-point properties for several 
parametrizations, characterized by the gauge parameters a, /3 
or 7 , as well as by the choice of the parametrization split pa¬ 
rameter r with r = 0 corresponding to the linear split ( 10 ) 
and r = 1, being the exponential split (11). Whether or not a 
field redefinition (30) is performed is labeled by “fr” or “nfr”, 
respectively. 


result would be exactly that of [62[ and in close agree¬ 
ment with [76, 77[ with slight differences arising from the 
regularization scheme. We summarize a selection of our 
quantitative results in Table I. 


B. Exponential split without field redefinition 

As a somewhat contrary example, let us now study the 
case of the exponential split ( 11 ) also without field redef¬ 
inition (30). Again, we find a non-Gaufiian fixed point. 
The corresponding estimates for the universal quantities 
at this fixed point in the Landau gauge limit a = 0 are 
displayed in Fig. 2. At first glance, the results seem simi¬ 
lar to the previous ones with a stationary point at 7 = 0 . 
However, the product 5 *A* shows a larger variation on 
the order of 5% and the critical exponent even varies by 
a factor of more than 40 in the range 7 S [—2,2]. We 
interpret the strong dependence on the gauge parameter 
7 as a clear signature that these estimates based on the 
exponential split without field redefinition should not be 
trusted. 

In fact, the real part of the critical exponents, Re0, 
have even changed sign compared to the previous case 
implying that the non-Gaufiian fixed point has turned 
UV repulsive. Similar observations have been made in 
[15] for the harmonic Feynman-type gauge a = 1 = /3 and 
an additional strong dependence on the regulator profile 
function Rk{x) has been found. We have verified that 
our results agree with those of [15] for the corresponding 
gauge choice. In summary, this parametrization serves as 
an example that non-perturbative estimates can depend 
strongly on the details of the parametrization (even for 
seemingly reasonable parametrizations) and the results 
can be misleading. The good news is that a study of the 
parametrization dependence can - and in this case does 
- reveal the insufficiency of the parametrization through 
its strong dependence on a gauge parameter. 
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Figure 2. Exponential split without field redefinition: resid¬ 
ual dependence of our estimates for the universal quantities 
on gauge parameter 7 in the limit a —>■ 0. A common sta¬ 
tionary point is again present at 7 = 0 , but the estimates 
for the universal quantities exhibit a substantial variation in 
the range 7 G [—2, 2]: g*A, varies by ~ 5% and Re 6 even by 
more than a factor of 40. The latter is a clear signal for the 
insufficiency of the parametrization. 


Figure 3. Linear split with field redefinition: residual de¬ 
pendence of our estimates for the universal quantities on the 
gauge parameter /? in the limit a —>■ 0. A common stationary 
point is approached for \j3\ — ^ 00 . Near the harmonic gauge 
/? = 1 (green dashed vertical line), both quantities have an 
extremum. For the whole range of /3 values, the estimates for 
the universal quantities exhibit rather small variations of 1 % 
for gr,A* and 10% for the more sensitive critical exponent Re0. 


C. Linear split with field redefinition 

For the remainder, we consider parametrizations of the 
fluctuation field which include held redefinitions (30). 
The canonical normalization achieved by these held re¬ 
definitions has not merely aesthetical reasons. An impor¬ 
tant aspect is that the nonlocal held redefinition helps to 
regularize the modes in a more symmetric fashion: the 
kinetic parts of the propagators then become linear in the 
Laplacian which are all equivalently treated by the regu¬ 
lator (31). A practical consequence is that the interplay 
of the degeneracy in the scalar sector no longer interferes 
with the regularization, i.e., the gauge parameter jS can 
now be chosen independently of a. Concentrating again 
on the Landau-gauge limit a —> 0, we observe for generic 
split parameter r that /3 = 0 no longer is an extremal 
point. 

Our estimates for the universal quantities for the case 
of the linear split (10) with field redefinition (30) and 
a —>■ 0 are plotted in Fig. 3. In order to stay away from 
the singularity at /? = 3, cf. Eq. (10), we consider values 


for /3 < 3 down to /I —>■ — 00 . As is obvious, e.g., from 
Eq. (21), the dependence of the propagator of the scalar 
modes and thus on /3 is such that the limits of large posi¬ 
tive or negative (3 —)■ ±c» yield identical results. Also the 
longitudinal ghost mode decouples in the limit /3 —?> ±00 
such that the whole flow in the large |/3|-limit is indepen¬ 
dent of the sign of /3. 

A non-Gaufiian fixed point exists, and a common ex¬ 
tremum of g*A* and Re0 occurs for /3 —>■ — 00 . Near 
(3 = 1 marking the harmonic gauge condition, both quan¬ 
tities are also close to an extremum (which does not oc¬ 
cur at exactly the same (3 value for both quantities). All 
fixed-point quantities for this case are listed in Tab. I (“fr 
T = a = 0, /3 = 1”). These values agree with the results 
of [37]. They are remarkably close, e.g., to those for the 
linear split without held redefinition. The situation is 
similar for the other extremum |,d| —> 00 (“fr r = a = 0, 
P = 00 ” in Tab. I). For the whole infinite /3 range stud¬ 
ied for this parametrization, (;*A* varies on the level of 
1%. The more sensitive critical exponent Re 6* varies by 
10% which is still surprisingly small given the simplic¬ 
ity of the approximation. Let us emphasize again that 
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Figure 4. Exponential split with field redefinition: residual 
dependence of our estimates for the universal quantities on 
the gauge parameter /3 in the limit a —>■ 0. A common sta¬ 
tionary point is approached for |/3| —^ oo, whereas no common 
minimum-sensitivity point is found near the harmonic gauge 
/3 = 1 or /3 = 0 (dashed vertical lines). Below j3 < —2, the 
critical exponents become real with the non-Gaufiian fixed 
point remaining UV attractive. For \j3\ —>■ oo, the results 
become independent of the gauge parameter a. 

varying (3 from infinity to zero corresponds to a com¬ 
plete exchange of the scalar modes from a (longitudinal 
vector component) to h (conformal mode) and hence to 
a rather different parametrization of the fluctuating de¬ 
grees of freedom. 

D. Exponential split with field redefinition 

Finally, we consider the exponential split (10), t = 1, 
with field redefinition (30). Having performed the latter 
has a strong influence on the stability of the estimates of 
the universal quantities at the non-Gaufiian fixed point, 
as is visible in Fig. 4. Contrary to the linear split, we 
do not find a common extremum near small values of (3: 
neither /3 = 0 nor the harmonic gauge (3 = 1 seem special, 
but, e.g., the product (;*A* undergoes a rapid variation 
in this regime. 

Rather, a common extremal point is found in the limit 
13 —>■ oo. In fact, becomes insensitive to the precise 
value of P for /3 < — 2 (with a local maximum near /? ~ 


—3, and an asymptotic value of g*A* ~ 1.033 for /3 —> 
oo. This estimate for 5 *A* is significantly larger than for 
the other parametrizations. The deviation may thus be 
interpreted as the possible level of accuracy that can be 
achieved in this simple Einstein-Hilbert truncation. 

As an interesting feature, the critical exponents be¬ 
come real for P < —2, and approach the asymptotic val¬ 
ues 9 = {4,2.148} for P —>■ 00 . The leading exponent 
9 = A reflects the power-counting dimension of the cos¬ 
mological term. This is a straightforward consequence 
of the fact that the A dependence in this parametriza¬ 
tion T = 1 , /3 —>■ 00 disappears from the propagators 
of the contributing modes. The leading nontrivial expo¬ 
nent 9 = 2.148 hence is associated with the scaling of 
the Newton constant near the fixed point, which is re¬ 
markably close to minus the power-counting dimension 
of the Newton coupling. The latter is a standard result 
for non-Gaufiian fixed points which are described by a 
quadratic fixed-point equation [87, 88 ]. The small differ¬ 
ence to the value 9 = 2 arises from the RG-improvement 
introduced by the anomalous dimension in the threshold 
functions (“ 77 -terms” as discussed in the Appendix). Ne¬ 
glecting these terms, the estimate of the leading critical 
exponents in dimension d is d and d — 2, as first dis¬ 
cussed in [17]. Also our other quantitative results for the 
fixed-point properties are in agreement with those of [17] 
within the same approximation. 

The significance of the results within this parametriza¬ 
tion is further underlined by the observation that the re¬ 
sults in the limit /3 —> 00 become completely independent 
of the gauge parameter a. In other words, the choice of 
the transverse traceless mode and the a mode (/3 —>■ 00 ) 
as a parametrization of the physical fluctuations removes 
any further gauge dependence. 

The present parametrization has also some relation to 
[89, 90], where in addition to the exponential split the 
parametrization was further refined to remove the gauge- 
parameter dependence completely on the semi-classical 
level. More specifically, the parametrization of the fluc¬ 
tuations was chosen so that only fluctuations contribute 
that also have an on-shell meaning. In essence, this re¬ 
moves any contribution from the scalar modes to the UV 
running. At the semi-classical level [89], the nontrivial 
critical exponent is 2 as in [17] and increases upon in¬ 
clusion of RG improvement as in the present work. The 
increase determined in [90] is larger than in the present 
parametrization and yields 0 ~ 3 which is remarkably 
close to results from simulations based on Regge calculus 
[91, 92] . 

The present parametrization with |/3| —> cx) is also 
loosely related to unimodular gravity, as the conformal 
mode is effectively removed from the fluctuation spec¬ 
trum. Still, differences to unimodular gravity remain in 
the gauge-fixing and ghost sector as unimodular gravity 
is only invariant under transversal diffeomorphisms. It 
is nevertheless interesting to observe that corresponding 
FRG calculations yield critical exponents of comparable 
size [93, 94] . 
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A similar conclusion holds for the more sensitive criti¬ 
cal exponents. Real and imaginary parts of the complex 
pair are shown in the lower panel of Fig. 5. Starting from 
larger values of a, it is interesting to observe that the 
imaginary part Im 0 decreases with decreasing a. This 
may be taken as an indication for a tendency towards 
purely real exponents; however, at about a = 1 this 
tendency is inverted and the exponents remain a com¬ 
plex pair in between Feynman gauge and Landau gauge 
within the present estimate. 

In summary, we observe no substantial difference be¬ 
tween the results in Feynman gauge a = 1 and those of 
Landau gauge a = 0 in any of the quantities of interest 
for the linear split and with field redefinition. Our results 
show an even milder dependence on the gauge parameter 
in comparison to the pioneering study of Ref. [72] , where 
the regulator was chosen such as to explicitly lift the de¬ 
generacy in the sector of scalar modes in the limit a —>■ 0. 
The present parametrization hence shows a remarkable 
degree of robustness against deformations away from the 
a-priori preferable Landau gauge. Hence, we conclude 
that the use of Feynman gauge is a legitimate option to 
reduce the complexity of computations. 


Figure 5. Linear split with field redefinition: dependence of 
estimates for the fixed-point values (upper panel) and the 
critical exponents (lower panel) on the gauge parameter a 
and harmonic gauge condition j3 = 1. No qualitative and 
only minor quantitative differences are found for the Feynman 
gauge a = 1 in comparison to the Landau gauge a — 0. 


In fact, the present parametrization allows for a closed 
form solution of the RG flow as will be presented in 
Sect. IV G. 


E. Landau vs. Feynman gauge 

Many of the pioneering computations in quantum grav¬ 
ity have been and still are performed within the harmonic 
gauge /? = 1 and with a = 1 corresponding to Feynman 
gauge. This is because this choice leads to a number of 
technical simplifications such as the direct diagonaliza- 
tion of the scalar modes as is visible from the off-diagonal 
terms in Eq. (19). Concentrating on the linear split with 
field redefinition, we study the a dependence for the har¬ 
monic gauge /? = 1 in the vicinity of the Landau and 
Feynman gauges. The results for the non-Gaufiian fixed 
point values are shown in the upper panel of Fig. 5. In 
essence, the fixed-point values show only a mild variation 
during the transition from the Landau gauge a = 0 to 
the Feynman gauge a = 1. In particular, the decrease 
of is slightly compensated for by a mild increase of 
A*. Effectively, the observed variation is only on a level 
which is quantitatively similar to other parametrization 
dependencies, cf. Table 1. 


F. Generalized parametrizations 

Having focused so far mainly on the gauge-parameter 
dependencies for fixed values of the split parameter r, 
we now explore the one-parameter family of parametriza¬ 
tions for general r. For this, we use the Landau gauge 
a = 0 and take the limit |/3| —>■ oo, where the fixed-point 
estimates of all parametrizations used so far showed a 
large degree of stability. Figure 6 exhibits the results for 
the non-Gaufiian fixed-point values (upper panel) and the 
corresponding critical exponents (lower panel). 

A comparison of the results for r = 0 and r = I reveals 
the differences already discussed above: an increase of 
the fixed-point values and the occurrence of real critical 
exponents for the exponential split r = 1. From the 
perspective of the principle of minimum sensitivity, it is 
interesting to observe that the fixed-point values develop 
extrema near r ~ 1.22. The product is maximal 
for T = 1 + ^ Also for this parametrization, 

the critical exponents of the fixed point are real and still 
close to the values for the exponential split, cf. Table 1. 
For even larger values of r, the critical exponents form 
complex pairs again. 

To summarize, in the full three-parameter space de¬ 
fined by r, /? and a > 0, we find a local extremum, 
i.e., a point of minimum sensitivity, at a = 0, /3 —)■ oo 
and r near the exponential split value r = 1. From this 
a-posteriori perspective, our results suggest that the ex¬ 
ponential split (with field redefinition) in the limit where 
the scalar sector is represented by the a mode may be 
viewed as a “best estimate” for the UV behavior of quan¬ 
tum Einstein gravity. Of course, due to the limitations 
imposed by the simplicity of our truncation, this conclu- 
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Figure 6. Parametrization dependence of fixed-point values 
(upper panel) and critical exponents (lower panel) as a func¬ 
tion of the split parameter r for the Landau gauge a = 0 
and I/?I —>■ oo. The fixed-point values exhibit extrema near 
T ~ 1.22, for the product of fixed-point values, this occurs at 
^ ~ ^ dashed vertical line). In this regime, 

the critical exponents are real and close to their values for the 
exponential split r = 1 (green dashed vertical line). 


sion should be taken with reservations. The resulting RG 
flow for r = 1 is in fact remarkably simple and will be 
discussed next. 


G. Analytical solution for the phase diagram 


Let us now analyze more explicitly the results for the 
RG flow for the exponential split with held redefinition 
in the Landau gauge and in the limit |/3| —>■ oo. Several 
simplifications arise in this case. The exponential split 
removes any dependence of the transverse traceless and 
vector components of the propagator on the cosmological 
constant. The remaining dependence on A in the confor¬ 
mal mode is Anally removed by the limit |/3| —> oo. As a 
consequence, the cosmological constant does not couple 
into the flows of the Newton coupling nor into any other 
higher-order coupling. Still, the cosmological constant is 
driven by graviton fluctuations. As emphasized above, 
any remaining gauge dependence on the gauge parame¬ 
ter a drops out of the flow equations. For the RG flow of 
Newton coupling and cosmological constant, we And the 


simple set of equations: 

43 810 \ 

47r 727r — 5g ) 

(35) 


A = = -2 - 


135g 
727r — 5g 


X-g 


In addition to the Gaufiian Axed point, these flow equa¬ 
tions support a Axed point at 


9* — 


1447r 
145 ’ 


A:,s 


48 

145’ 


g^ A:,s 


691271 

21025’ 


(36) 


cf. Table 1. Also the critical exponents 6i being 
(—1) times the eigenvalues of the stability matrix 
9P{g,\)/d{g,X) can be determined analytically. 


do = 4, 01 


58 


(37) 


The fact that the largest critical exponent corresponds to 
the power-counting canonical dimension of the cosmologi¬ 
cal term is a straightforward consequence of the structure 
of the flow equations within this parametrization: as we 
have g = {2 + g{g))g and A = (—2-|-?7(g))A-|-(!l(g), the ex¬ 
istence of a non-Gaufiian Axed point requires g(g*) = — 2. 
As the stability matrix is triangular, the eigenvalue asso¬ 
ciated with the cosmological term must be —4 and thus 
Oq = 4. Rather generically, other parametrizations lead 
to a dependence of rj also on A and thus to a more in¬ 
volved stability matrix. 

In the physically relevant domain of positive gravita¬ 
tional coupling g > 0, the Axed point g* separates a 
“weak” coupling phase with g < g* from a “strong” cou¬ 
pling phase g > g*. Only the former allows for trajec¬ 
tories that can be interconnected with a classical regime 
where the dimensionless g and A scale classically, i.e., 
g ~ 2g and A ~ — 2A such that their dimensionful coun¬ 
terparts approach their observed values. Trajectories in 
the strong-coupling phase run to larger values of g and 
terminate in a singularity of fig at gsing = 727r/5 indicat¬ 
ing the break-down of the truncation. 

All trajectories in the weak coupling phase with g < g* 
run towards the Gaufiian Axed point for g and thus, also 
the flow of A in the infrared is dominated by the Gaufiian 
fixed point. This implies that all trajectories emanating 
from the non-Gaufiian fixed point with g < g* can be 
continued to arbitrarily low scales, i.e., are infrared com¬ 
plete. They can thus be labeled by their deep infrared 
value of gA approaching a constant, which may be iden¬ 
tified with the product of Newton coupling and cosmo¬ 
logical constant as observed at present. A plot of the re¬ 
sulting RG flow in the plane (g, gA) is shown in Fig. 7. It 
represents a global phase diagram of quantum gravity as 
obtained in the present truncation/parametrization. We 
emphasize that no singularities appear towards the IR 
contrary to conventional single-metric calculations based 
on the linear split. 
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Figure 7. Global phase diagram in the {g,gX) plane for the exponential split with field redefinition and |/3| —>■ oo. Arrows 
point from IR to UV indicating the approach to the UV fixed point at g* = 1447r/145 and A* = 48/145. The color indicates a 
measure for the flow velocity, (dtg)^ + {dt{gX /\/l + g^X^))^. 


The flows (34) and (35) can be integrated analytically. 
Converting back to dimensionful couplings, the flow of 
the running Newton coupling G{k) satisfies the implicit 
equation, 


Gn = 


G{k) 


(l-^PG(fc)) 


27 ? 
29 


(38) 


where Gw is the Newton coupling measured in the deep 
infrared fc —> 0. Expanding the solution at low scales 
about the Newton coupling yields 

G(fc)~Gw(^l-^fc2Gw + G((A:2Gw)")^ (39) 

exhibiting the anti-screening property of gravity. 

The flow of the dimensionful running cosmological con¬ 
stant A(fc) can be given explicitly in terms of that of the 
running Newton coupling, 

A{k) = -h (I447r - 145G(fc)fc2) ^ 

25 lOTT 

1 4477 / ^ \ 

( 40 ) 


Here, i = (125AGw + 432^), 

and A is the value of the classical cosmological constant 
in the deep infrared fc —>■ 0. The low-scale expansion 
about fc = 0 yields 

A(fc) ~ A (^1 - ^k^GN + O (^^AGw, (A^^Gw)')) 

(41) 

Thus, A{k)/G{k) = A/Gn + 0{k^), implying a compar¬ 
atively slow running of the ratio towards the UV. This 
explicit solution of the RG flow might be useful for an 
analysis of “RG-improved” cosmologies along the lines of 
[95-101]. 


H. Generalized ultra-local parametrizations 

For the most general, ultra-local parametrization (13), 
it turns out that the flow equation in our truncation does 
only depend on the linear combinations Ti := r/4 -|- ra 
and T 2 := 72/4 -1- T 4 , leaving only two independent 
split parameters. Instead of exploring the full high¬ 
dimensional parameter space, we try to identify relevant 
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points as inspired by our preceding results. For instance 
for the choice Ti = 1/4 and T 2 = —1/8, any dependence 
on a drops out, indicating an enhanced insensitivity to 
the gauge choice. The resulting flow equations are 


5 = 23 + 


135(^ - 3)32 


A = -2A + 


(5/3 - 3)3 - 72(/3 - 3)^ ’ 

3((-669 + 215/3)3 + 36(/3 - 3)7r(4 - 15A)) 


(42) 


47r((3-5/3)3 + 72(/3-3)^) 


(43) 


In the limit \P\ —)■ 00, these are identical to the exponen¬ 
tial split in the same limit. The non-Gaufiian fixed point 
occurs at 


1447r(/3 - 3) ^ 48(/3 - 3) 

“ 145/3-411’ “ 145/3-411’ 

3 *A* = 69127r \ •—— -• | . (44) 

^ V145/3-411 ) ^ ’ 

Apart from the pathological choice /3sing = 3 (incom¬ 
plete gauge fixing) where this fixed point merges with 
the Gaufiian fixed point, no further extremal point is ob¬ 
served except for the limit |/3| —>■ 00 . The critical expo¬ 
nents are 


^0 — 4, 


58 16 

^ 45(^ - 3) ■ 


(45) 


Also the exponents become minimally sensitive to the 
choice of /3 for |/3| —> 00 . 

As an oddity, we mention the particular case /3 = 3/5, 
where the flow equations acquire a pure one-loop form. 
In this case, the second critical exponent is exactly 2 as 
it must. 

More importantly, the interdependence of gauge and 
parametrization choices is also visible in the following 
fact: we observe that the choice of the gauge parameter 
1/3 1 —> cx) removes any dependence of our flow on the 
parameter T 2 independently of the value of a. In other 
words, this limit brings us back exactly to the case which 
we discussed above in Sect. IVF, such that the seemingly 
much larger class of parametrizations (13) collapses to a 
one-parameter family. 


I. Arbitrary dimensions 

Finally, we discuss the stability of the UV fixed-point 
scenario and its parametrization dependence in arbitrary 
dimensions, focusing on d > 2 (for a discussion of d = 2 
in the present context, see [15, 17, 89]). In fact, there are 
some indications in the literature that the parametriza¬ 
tion dependence is pronounced in higher dimensions. 
Whereas standard calculations based on the linear split 
generically find a UV fixed point in any dimension d > 2 
and gauge-fixing parameter a, see e.g. [102, 103[, a recent 
refined choice of the parametrization to remove gauge- 
parameter dependence on the semi-classical level arrives 


at a different result [89, 90[: the UV fixed point can be 
removed from the physical region if the number of phys¬ 
ical gravity degrees of freedom becomes too large. As 
the latter increases with the dimensionality, there is a 
critical value dcr above which asymptotically safe grav¬ 
ity does not exist. The resulting scenario is in line with 
the picture of paramagnetic dominance [104, 105[, which 
is also at work for the QED and QGD /3 functions: the 
dominant sign of the /3 function coefficient arises from the 
paramagnetic terms in the Hessian which can be reversed 
if too many diamagnetically coupled degrees of freedom 
contribute. 

Our results extend straightforwardly to arbitrary di¬ 
mensions. Starting, for instance, with the most general 
parametrization (13) in d dimensions, the flows of 3 and 
A depend only on the linear combinations Ti = r/d -I- T 3 
and T 2 = T 2 /d + T^. Gomparable results as in d = 4 di¬ 
mensions apply: in the limit of |/3| —> 00, also T 2 drops 
out such that a one-parameter family remains. In turn, a 
complete independence of the gauge parameter a can be 
realized with the parametrization specified by Ti = 1 /d 
and T 2 = -l/( 2 d). 

We illustrate the stability properties of the asymptotic- 
safety scenario in arbitrary dimensions by choosing the 
Landau-gauge limit a —>■ 0 as well as |/3| —>■ 00 , keep¬ 
ing Ti as a free parameter. Then, we know a priori that 
Ti = 1 /d would be a preferred choice from the view point 
of gauge invariance; it would also correspond to the expo¬ 
nential parametrization r = 1, ra = 0. Fig. 8 displays the 
fixed-point values for 3 * A* as a function of Ti for various 
dimensions d = 3,..., 7. While d = 3 exhibits a rather 
small parametrization dependence, d = 4 reproduces the 
earlier results of Fig. 6 (upper panel) now as a function 
of Ti with an extremum not far above Ti = 1/4. By 
contrast, 3 *A* develops a kink for d = 5 that turns into 
a singularity for d = 6 and larger. For increasing d, the 
kink approaches the preferred parametrization Ti = 1 /d 
(vertical dashed lines in Fig. 8 ). The singularity in 3 *A* 
occurs for a critical dimension dcr — 5.731. 

This observation suggests the following interpretation: 
whereas we can identify a UV fixed point for any dimen¬ 
sion as long as we choose Ti sufficiently far away from 
Ti = 1 /d, we find a stable fixed-point scenario only for 
d = 3 and d = 4 integer dimensions. Already for d = 5, 
the fixed-point product 3 * A* can change by two orders of 
magnitude by varying the parametrization, which is at 
least a signature for the insufficiency of the truncation. 
For d > dcr — 5.731, 3 *A* can become unboundedly large 
as a function of the parametrization, signaling the insta¬ 
bility of the fixed point. 

If these features persist also beyond our truncation, 
they suggest that the asymptotic safety scenario may 
not exist far beyond the spacetime dimension d = 4. 
Whereas this does not offer a dynamical explanation of 
our spacetime dimension, it may serve to rule out the mu¬ 
tual co-existence of extra dimensions and asymptotically 
safe quantum gravity. 








14 


10 ®. 


1000 



0.1 


-0.1 0.0 0.1 0.2 0.3 0.4 0.5 


Figure 8. Parametrization dependence of fixed-point value 
for g*A* as a function of the split parameter T\ in the Lan¬ 
dau gauge Of = 0 and |/3| —>■ oo for different dimensions 
d = 3,4, 5, 6, 7 (from bottom to top). Vertical lines mark the 
value of the parameter T\ = Ijd preferred by gauge-parameter 
a independence. For d > dcr — 5.731, the fixed-point product 
g*A* develops a singularity at Ti = 1/d. 


terplay between various aspects of parametrizations, e.g., 
gauge choice vs. field parametrization, can arise. With 
hindsight, the results obtained within the exponential 
split with field redefinition in the limit where the graviton 
degrees of freedom are spanned by the transverse trace¬ 
less and a scalar mode (|/3| —>■ oo) exhibit the highest 
degree of comprehensiveness: complete independence of 
the gauge parameter a, fully analytical and integrable 
global flows with a classical IR limit in the physical pa¬ 
rameter regime, real critical exponents at the UV fixed 
point, and the existence of an upper critical dimension 
for the asymptotic safety scenario. The exploration of 
higher-order truncations [106] and the inclusion of mat¬ 
ter degrees of freedom [107-109] in this parametrization 
appears highly worthwhile, c.f. [17, 18] for scalar matter. 

In summary, our work exemplifies that a careful in¬ 
vestigation of parametrization dependencies facilitates 
both a test of the robustness of nonperturbative quan¬ 
tum gravity computations as well as the identification of 
a parametrizations which may be better adapted to the 
physical mechanisms. 


V. CONCLUSIONS 
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Our results show a remarkable stability in a variety of 
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rather large stationary regimes in the space of all param¬ 
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two RG relevant directions, corresponding to the New¬ 
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cal parameters. For most parametrizations, the universal 
quantities show a remarkably mild (given the simplicity 
of the approximation) variation and thus a high degree of 
stability. Our scan of parametrization dependencies can 
also help identifying less robust parametrizations, and 
thus help judging the physical relevance of results. 
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Appendix A: Flow equations 


In this section, we display the right hand side of the 
Wetterich equation for general Regulators i?fc[A] and in 
dimension d = 4. For simplicity, we introduce the anoma¬ 
lous dimension Tj = [g — 2g) /g, and refer to terms linear 
in g as “r^-terms”. Let us start with the contribution from 
the TT-mode: 


5 ^^ = 2^2 




-12^ 


+ {4 - 3t)Q2 


A-k R/c - 2A(1 - r) 


A + Rfe-2A(l-r) 

Rfc V^k 

(A + Rfe-2A(l-r))2 


(Al) 


where the Q functionals are defined in terms of Mellin 
transforms [20] . For the transverse vector, and without 
field redefinition, let us define 


riT _ 
= 


-{Rk - gRk) (^2A(1 - r) - ^(Rfe + 2A)^ 
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-2(A + 2A)(l-r)i?fe X 


(A + Rk) ( ^^+2A(1 -t) 


(A2) 


With that, we have 


For the scalar contribution, we first define 


TT-" = -— (4aA(l - t)(A + Rk)^ + (a _ 3)(A + Rkf) 


= ^Q2 [^r] + R Qqi [^r] (A3) 

+ 3(1 _ ^(1 _ _ 3^(1 _ _ 


On the other hand, the contribution with field redefini¬ 
tion reads. 


tt'' = - (- 4aA(l + t) -h (3a - /3^)(A -t Rk)) (A6) 

7T- = -^ia-P){A + Rkr 


(A5) 

(A6) 

(A7) 


P = 


4a 


(4a(A + (2 - 77)A)(1 - t)( 2A + Rk)Rk 


= ^Q2 


Rk — rjRk 


A + Rk — 2aA(l — t) 


+ 8aA(l — t)(A 4- Rk)Rk 

+ 3(a — 3)(A -|- Rk)^Rk 

-(a - 3)(3A2 + iARk + RDpRk) 

p'^ = -T^(3a - /3^)(7?fe - pRk) 
loa 


(A8) 

(A9) 


i'Q. 

Rk — pRk 


A + Rk — 2aA(l — r) 


(A4) 


P" = -^(a - ^)(2i?fc(A + Rk) - pRki 2 A + Rk)). 

8a 


(AlO) 


-|--(1 - 2a(l - t))Q 2 


Rk pRk 


{A + Rk - 2 aX{l - t))^ 


The contribution is 


5"" =^Q2 


TT'^p*' -h - 2tT^P^ 




+ R< ^Qi 




- ^(6- a(4- 3 t ))(34 


t(T /Trh 




+A(1 — t)Q 3 


+ —(6 - a(4 - 3 t))Q 4 
4a 


P^ 

--0, 

P"" 


p^ 


32^" 


4a 



- 2Tr^p^) 


T 

- Q 2 

32^^ 


^TT'^TT^ — 

-h - 2tt^p^) 


+ -^(1 ~ "T)Qz 


TT*' -t - 27r’^p’^) 


^7r°'7r*' — 


a — j: 
16a 


Q 2 


(TT-^ph -t ttI'p'^ - 27r’^p’^) 


With field redefinition, define 


TTfr = -— (4aA(l -t) + (a - 3)(A -I- Rk)) 

4a 

TT^ = ~Y5^ (—4aA(l -f r) -l- (3a — 0 ^){A -\- Rk)) 

nl = -^ia-(3)i^ + Rk) 

Pfr = ^(3 - a)(i?fc - pRk) 

Pfr = “Y^(3a - l 3 ‘^)iRk - pRk) 

Pfr = - pRk). 


(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 


Then, the scalar contribution is 


1 


= 2^2 


7rgpg.+ ttIpI - 27rj(.p| 

- iKf 


R\^Qi 


TTfrPfr + TTfrPfr " 27r|p| 
TT^TT^ - (TTg.)^ 


- ^(l-a(l-T))Q 2 


Pfr 


^TTfr - (TTg.) 
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T 

— Q 2 

32^ 


T 

-<32 

32^ 


Ph 


a —13 
16a 


Q 2 


KK - (^fr) 

- (’^fr)") 


Pi, 


kk - 

a —13 


+ ^(l-a(l-T))(32 


’"fr {KPh + ^hPi, - ‘^'^IpI) 


(^TT^ - (TTg.)^) 


16a 


<32 


{KPh + ^frPfr - ‘^^IpI. 


(A18) 


Further, the ghost contribution reads without field redef¬ 
inition, 




Rk 


(A + i?fc)2 


(A20) 


= - 5 Q 2 


+j«^ 


With field redefinition, it is 

Rk 


Rk 

-«(> 

Rk 

A J- Rk 

A J- Rk 

Rk 

1 0, 

Rk 

(A + 


(A + Rk) 


(A19) 


= - 4 Q 2 


A -|- Rj. 




Rk 


A + 


Finally, the contribution of the Jacobian for the case 
without field redefinition is 


^Jac__ ^gh I 


?=0 


Q 2 


Rk 


A + Rk 


,RQi 


Rk 


(A + i?fe)2 

(A21) 


With field redefinition, all Jacobians cancel, at least on 
maximally symmetric backgrounds, which is sufficient for 
the truncation considered here [72]. 
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